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Abstract. - Formation of turbulence of capillary waves is studied in laboratory experiments. The 
spectra show multiple exponentially decreasing harmonics of the parametrically excited wave which 
nonlinearly broaden with the increase in forcing. Spectral broadening leads to the development 
of the spectral continuum which scales as oc /~ 2 ' 8 , in agreement with the weak turbulence theory 
(WTT) prediction. Modulation instability of capillary waves is shown to be responsible for the 
transition from discrete to broadband spectrum. The instability leads to spectral broadening of the 
harmonics, randomization of their phases, it isolates the wave field from the wall, eventually allows 
the transition from 4- to 3-wave interactions as the dominant nonlinear process, thus creating the 
prerequisites assumed in WTT. 



Broad frequency and wave number spectra character- 

• ^nize many turbulent systems. When energy is injected into 
^a system in a limited range of wave numbers or frequen- 
(— | cies, nonlinear interactions extend this range. In the flows 
O described by the Navier-Stokes equations such spectral en- 

1 ergy transfer at large Reynolds numbers is attributed to 

^^the energy cascades. The idea of the cascade also lies in 
£>the center of the weak wave turbulence theory (WTT) 
^\| [l] , where instead of strong nonlinearity of the Navier- 

Stokes equation, the cascade is driven by weakly interact- 
^2 m § waves having almost random phases. Among the best 

known applications of the WTT to the surface waves is 
^Othe theory of capillary wave turbulence of Zakharov and 
^Filonenko |]. 

^""" l The theory assumes that the dispersion relation of capil- 
• «lary waves Wk = (ak 3 / p) 1 ^ 2 (where a is the surface tension 

• £h and p is the liquid density) allows simultaneous satisfac- 

tion of the matching rules for the wave triads in the fre- 
S-h quency and in the wave number domains " for any possible 



choice of the wave vectors" ki,k 2 ,k 3 



ki =k 2 + k 3 . 



(1) 
(2) 



Random-phase three-wave interactions are then consid- 
ered the main mechanism of the spectral energy transfer 
in capillary wave turbulence from low to high k and /. 
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The theory [2] predicts the wave number spectra of the 
surface elevation oc fc -19 / 4 , while the frequency spec- 
trum is given by 



E u 



F 1 /2 /9 -2/3 a l/6 w - 



17/6 



(3) 



The WTT has been tested in several experiments |4f|10| 
which gave contradicting results. Some papers report good 
agreement of the spectrum shape with theory [4]-[7] , while 
some others report deviations from the expected / or k 
power-law scalings [8 -10 



Since in many laboratory experiments surface waves are 
excited parametrically, it is important to clarify whether 
the assumptions of the WTT are satisfied there. In par- 
ticular, the assumptions are: (a) randomness of the waves 
phases, (b) infinite-box size, (c) the dominance of the 
three-wave interactions as the main nonlinear wave pro- 
cess. It is important to recall that the frequency spectra 
observed in experiments with the monochromatic para- 
metric forcing of capillary waves are often dominated by 
multiple harmonics of the first subharmonic of the forc- 
ing frequency (e.g. 
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These harmonics are ob- 
served starting from the threshold of the parametric ex- 
citation. At higher drive the spectra may broaden show- 
ing the power-law continuum (4]-[6 11 . It is not under- 



stood (a) how the harmonics are generated and (b) how 
the power-law continuum is formed. Understanding the 
above effects and assumptions is necessary for evaluation 
of (dis)agreement between experiments and theory. 
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Another important question which needs to be ad- 
dressed with regard to turbulence formation in parametri- 
cally driven capillary waves is what type of the wave-wave 
interactions is dominant in each of the regimes (low/high 
excitation drive). It is usually assumed that the three- 
wave interactions is the only dominant process here. Para- 
metric excitation can be viewed as a degenerate process 
of a three-wave decay of the infinitely long wave, which 
has a wave number kg = and a frequency lo , into two 
oppositely propagating subharmonic waves: 



k = ki + k 2 « 0, 



(4) 
(5) 



where k\ — —k^ and u)\ = ui2 — {ak\/ p) 1 / 2 . This leads to 
the generation of the waves at the frequency u)\ — w /2 
(Faraday waves). 

Apart from the decay of the fco = wave, simultaneous 
satisfaction of the Eqs.Q and Q is severely restricted for 



capillary (and capillary-gravity) waves, see e.g. |12||13 
The selection rules favor interacting triads with disparate 
wave vectors, for example, those with k\ « k 2 ,k 3 , and 
the angles between interacting wave vectors in a narrow 
band between 75° and 82°. In the presence of discrete 
spectra which contain only narrow multiple frequency har- 
monics, as is commonly observed, the above restriction 
practically forbids three-wave interactions unless spectral 
lines sufficiently broaden. However, before the broad spec- 
tra are formed, no three-wave processes are allowed, ex- 
cept for the parametric decay of the forcing. 

In this Letter we summarize recent findings and present 
new results on the formation of continuous frequency spec- 
tra of capillary waves during parametric excitation. At 
modest damping (in distilled water) parametrically ex- 
cited waves are unstable to small perturbations of the wave 
amplitude. This modulation instability seems to be re- 
sponsible for the randomization of the wave phases, the 
detachment of the wave field from the container walls, 
spectral broadening, and gradual transition from discrete 
to continuous spectra. 

The experimental setup is similar to that described in 
11 14 . Waves are excited parametrically in a vertically 



shaken container. Circular containers of the inner diame- 
ters > 15 cm filled with distilled water to the heights of > 
20 mm are mounted on the top plate of the 4k N electro- 
dynamic shaker. The container is shaken at the frequency 
of fo =(0-3500) Hz. The acceleration a can be varied in a 
broad range. The amplitude of the shaker vibrations Ao 
is proportional to the acceleration at a given oscillation 
frequency, Ao cx a. Below we will refer to the acceleration 
above the threshold of parametric excitation Aa = a — ath, 
rather than to its absolute value. 

The surface elevation of the fluid is measured using both 
the light diffusion and the laser reflection techniques de- 
scribed in 11 14 . In the light diffusion technique, skim 
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milk is added to the water as a diffusing agent as proposed 



Fig. 1: Snapshots of the wave field evolution during the 
startup: (a) t=0.125 s, (b) £=0.25 s, (c) t=0.625 s. Autocor- 
relation functions of the surface elevation during parametric 
excitation at fo = 60 Hz at the accelerations of (d) Aa = 0.25p, 
(e) Aa = 0.6g , (f) Aa = 1.3s . 



in A thin laser beam is shined vertically from below 
the cell. A diffusive screen is mounted above and near the 
surface to collect the transmitted light (whose intensity is 
proportional to the surface height), which is imaged into 
a photo-multiplier tube. The transmission of a thin laser 
beam reduces the problem of the surface averaging (the 
"cross-over length" effect discussed in In our setup 
the diameter of the probing laser beam at the fluid sur- 
face (d,B ~ 0.2 mm ) is substantially smaller than typical 
wavelengths of capillary waves in the frequency range be- 
low 2 kHz (A ~ 0.5mm) such that the measured spectra 
are not affected by the spatial averaging and the inten- 
sity of the transmitted light remains proportional to the 
wave height rj. As a result, measurements can be com- 
pared directly with the theory predictions. At low wave 
amplitudes (e.g. at the high frequency excitation) we use a 
more sensitive technique based on the reflection of a thin 
laser beam (< 0.5 mm). This technique is not used at 
high forcing levels to avoid signal clipping by steep waves. 
Spectra obtained using data from both techniques are very 
similar. 

First, we consider the dynamics of the wave excitation 
captured using fast video camera. Figures [lja-c) show 
three snapshots of the wave field filmed at the intermediate 
acceleration in a cylindrical container shaken vertically at 
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the frequency of 250 Hz. Initially, 0.125 s after the forcing 
is switched on, Fig. [lja), the wave fronts of the excited 
wave are concentric circles, reflecting the shape of the wall. 
However very soon these wave fronts become azimuthally 
modulated, Fig.JlJb), and after about 0.6 s the wave field 
appears to be severely disordered, as in Fig. We). Similar 
evolutions are observed for all excitation frequencies. 

Order-disorder transitions in capillary waves have been 
studied in 



E(f) (a.u.)l 
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where it has been shown than the radial 
correlation length is dramatically reduced at rather low 
excitation levels. When the transition occurs, the bulk 
wave field is no longer connected to the wall since the cor- 
relation length becomes much shorter than the container 
radius. At this point hnite-box size effects become unim- 
portant. Similarly, we observe a dramatic change in the 
autocorrelation time with the increase in the acceleration. 
Figures |T|(d-f ) show the autocorrelation functions of the 
surface elevation measured at the mid-radius at three lev- 
els of accelerations. At low acceleration, the autocorrela- 
tion function shows two characteristic frequencies: 30Hz 
of the first subharmonic of the excitation frequency, and 
a modulation at about 4-5 Hz corresponding to the eigen- 
mode of the cylindrical container (wave length of 14 cm, 
which is the container radius), Fig. [ljd). At higher forc- 
ing, the autocorrelation function becomes exponentially 
decaying, as in Fig. |T](e). Further increase in the acceler- 
ation leads to a fast decrease in the auto correlation time 
down to 1-2 periods of the first subharmonic wave (30 Hz) 
illustrated in Fig. [ijf). 

Such a transition from order to disorder is due to the de- 
velopment of the modulation instability of capillary waves. 
This instability, which was first considered theoretically 
by Zakharov in 1968 16 , is the capillary counterpart of 



the well-known Bcnjamin-Fcir instability of gravity waves 



17 . Modulation instability (MI) of capillary waves has 
been identified for the first time in 111]. It develops when 
the Lighthill criterion, (duo / d\a\ 2 ){(TUJk/ dk 2 ) < 0, (where 
lo = Wfe[l — (fca) 2 /16] includes a nonlinear frequency cor- 
rection [18]) is satisfied. The development of the MI is 
observed at a given point in space as the modulation of 
the envelope which, at higher wave amplitudes (forcing) 
leads to breaking of the initially continuous waves into 



sequence of the envelope solitons 11 . 

To understand how and when the transition to the 
power-law continuum occurs, we now consider frequency 
spectra measured at various levels of acceleration. Fig- 
ure [2] shows frequency spectra measured at three different 
accelerations, up to the threshold of the droplet gener- 
ation. At the lowest acceleration, just above the para- 
metric excitation threshold, we observe several (up to 10) 
frequency harmonics f n — nfi which become broader 
as the acceleration is increased. To obtain a spectral 
power contained in each of the harmonics, spectral den- 
sities E(f) are integrated within the frequency bands 
[fn — /i/2, f n + fi/2}. The resulting spectral powers Ef are 
shown by open diamonds above the corresponding spectral 
peaks in Figs. [2]ja-c). 
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Fig. 2: Frequency spectra of capillary waves at different ac- 
celerations, (a) Aa = 0.5g, (b) Aa — lAg and (c) Aa — 2.1g. 
Open diamonds show the spectral powers Ef of each harmon- 
ics. 



At low forcing (Aa = 0.5g), Fig.[2ja), the harmonic am- 
plitudes exponentially decay as a function of the wave har- 
monic, E(nfi) ~ e~ n . The band- integrated amplitudes of 
the harmonics in Fig. [2^a) also follow the geometric pro- 
gression law Ef — E\e~~ A f I f 1 , where A — 1.2. Along with 
the harmonics, a zero-frequency sideband (ZFSB) is ob- 
served at all levels of forcing, even at the lowest one (/ < 1 
Hz). Both the harmonics and the ZFSB are observed just 
above the threshold of the parametric excitation of the 
first harmonic, at Aa > 0. The ZFSB results from the 
modulation instability of capillary waves. 

The generation of the multiple frequency harmonics, 
however is not really understood. Linear analysis of the 



stability of the interface between two fluids 19 shows that 
the parametric excitation leads to the instability zones 
corresponding to the multiple frequency harmonics and, 
in principle, could explain the harmonic generation. How- 
ever, the acceleration thresholds for the generation of har- 
monics with n > 4 expected for water are much higher 
(20 — 30<7 a t fo = 50Hz) [20] than observed in our ex- 
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Table 1: Energy distribution in various spectral bands cor- 
responding to spectra of Fig. [2j^ Energy is shown as a frac- 
tion of the total energy Et = J z E(f)df. Ezf, Eex P and 
Ep ow refer to the zero-frequency sideband, exponential, and 
the power-law ranges respectively. 



periments (0.2 — 0.5g at f = 60Hz). Another possible 
interpretation is related to the WTT energy flux, which, 
as was shown in numerical simulations 



21 can lead to 



the accumulation of energy in the harmonics of the forc- 
ing frequency. This however is unlikely to be the case in 
this experiment since the harmonics are observed at very 
low forcing levels, when all wave harmonics are perfectly 
coherent 
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which is inconsistent with the assumption 
of the random phases. 

An intermediate regime is achieved at the forcing level 
of Aa = lAg. Fig. [2jb) shows the spectrum in this case. 
The spectrum is also dominated by discrete harmonics, 
however spectral broadening of the harmonics leads to 
their overlapping and to the formation of spectral con- 
tinuum. It can be seen that up to the frequencies of about 
220 Hz, the spectral power Ef follows an exponential fit, 
Ef = E\e~ A 5l$ x (A = 0.9), while at higher frequencies 
the spectrum scales according to a power law, Ef ~ /~ 3 ' 2 . 
The ZFSB in this regime extends to about 10 Hz. 

At the highest level of forcing (Aa = 2.1g), Fig. [2jc), 
the spectrum changes further, showing exponential range 
Ef — E%e~ A f l^ 1 (A = 0.8) for the first 4 harmonics and an 
extended continuum which is well fitted by a f~ 2 8 power 
law for the frequencies above 120 Hz. The E(f) ~ /~ 2 ' 8 
part of the spectrum, which extends over one decade, is in 
a good agreement with the WTT E ~ w~ 17 / 6 prediction, 
Eq. ^. The energy contained in the ZFSB in this regime 
is the highest. 

Table 1 shows the energy distribution between three 
spectral bands identified above: ZFSB, exponential (dis- 
crete harmonics) and the power-law (continuum) ranges. 
The forcing injects energy into the exponential (discrete 
harmonic) spectral range, which is then redistributed into 
other ranges via wave- wave interactions. At low levels of 
forcing (Aa = 0.5g) no transfer of energy into the high- 
frequency continuum is observed, while 18% of the total 
spectral energy is coupled into the ZFSB. At the interme- 
diate forcing, corresponding to Fig.[2jb), only 1% of total 
energy is contained in the power-law continuum, while the 
rest is distributed almost equally between ZFSB and the 
exponential range. Finally, at the highest forcing, 23% of 
energy is delivered from the forcing range into the power- 
law continuum, leaving 31% in the ZFSB. 

This can be explained by recalling that at low forcing 
no three-wave interactions between discrete harmonics are 



possible. As the drive is increased, the satisfaction of the 
three- wave matching rules of Eqs.([T]) and ^ in this sit- 
uation becomes possible due to the spectral broadening 
of the harmonics, when broadening compensates for the 
frequency mismatch Afl. This broadening is due to the 
modulation instability, a four-wave process of the side- 
band generation. Thus, a four-wave process leads to the 
nonlinear broadening of the harmonics and to the genera- 
tion of the ZFSB until broadening exceeds AO such that 
three-wave interactions arc allowed. This opens the gate 
for the energy flux into the power-law continuum. The 
power law tail Ef ~ /~ 2 ' 8 , observed at the highest forc- 
ing, is in agreement with the the WTT scaling of w -17 / 6 . 
The random-phase assumption is justified at high forcing 
as was reported in 
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Also, the build-up of the spectral 
continuum increases the range of waves satisfying Eqs.Q 
and ^ and opens the way for the three- wave interactions. 

The evidence in support of the four-wave interaction 
process in the harmonic-dominated regime has been pre- 
sented in 
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If the wave harmonics are modulation- 
ally unstable, each of them satisfies four-wave synchro- 
nism conditions: w„_ + (J n + = 2ui n , k\ + &2 = 2k n , or in 
another form, w n _ -Tl+ = uj n ± il, k\ y 2 = k n ± K, where 
f2 is the modulation frequency 22 . In this case the wave 
harmonics u) n are all phase coupled via Q, represented in 
Fig. [2] by the ZFSB range. The coherent-phase four-wave 
coupling has indeed been found in this system. As has 
been reported in [14], the maximum tricoherence, which 
is a degree of a coherent four-wave coupling, is high in the 
conditions corresponding to the discrete spectra of capil- 
lary waves. This result also indicates that the four-wave 
process (cubic nonlincarity) is at work in the conditions 
when the three-wave interactions (quadratic nonlinearity) 
is forbidden by the Eqs.Q and (2}. 

Let us now analyze physical conditions which determine 
the transition of the exponential part of the spectrum into 
the power-law continuum. We denote the cross-over fre- 
quency where these two ranges meet as f co . It can be seen 
from Figs. [2jb,c) that f co decreases with the increase in 
forcing from f co w 220 Hz at Aa = lAg to f co w 120 Hz 
at Aa = 2.1g. We have previously reported that the spec- 
tral line shapes of the broadened wave harmonics can be 
accurately represented using the hyperbolic secant 11 : 
E n oc sech[b n (f — f n )]. In the time domain such spec- 
tra correspond to the strongly modulated sinusoidal pulses 
s n (t) = (7r /b n )sech[TT 2 1 (b n t)]e l f nt which are very similar 
to the well known envelope solitons. Statistically aver- 
aged lifetime r„ (length) of such a pulse is related to the 
spectral width parameter b n as r„ = b n /iT 2 . The loss 
of coherence in the wave-wave interactions occurs when 
the lifetime r„ of a highest-frequency interacting wave be- 
comes shorter than the period of the first harmonic T\\ 
& t T n < T\. We find that exactly this condition is sat- 
isfied at / = f co : t co R3 T\. Figure [3] shows time evo- 
lution of the two wave harmonics bandpass-filtered from 
the frequency bands [/„ - fi/2, /„ + fi/2], for /i = 30 Hz 
(n = 1) and /g = 240 Hz (n = 8) for the data correspond- 
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Fig. 3: Time evolution of the n = 1 frequency harmonic, / = 
(30 ± l5)Hz, (top) and of the n = 8 harmonic at / = (240 ± 
15)Hz for the intermediate forcing regime of Fig. [2jb) . 



ing to Fig.|2|b). The second frequency is chosen such that 
/s ^ fco ~ 220 Hz. One can see that the average length 
of the /g pulse, r, is indeed comparable to the period of 
the first harmonic T\ = 1/fi- This result is confirmed 
by calculating b n from the spectra at f n > f co at several 
levels of forcing and by comparing t„ = b n /ir 2 with T\. 

Naturally, when r n < T± , the nonlinear interaction time 
between the 1st and the n-th harmonics decreases dramat- 
ically turning coherent interaction into the random-phase 
one. Thus, we conclude here that the criterion correspond- 
ing to the transition from the exponential to the power-law 
range coincides with the criterion of phase-randomization 
in the wave- wave interactions, namely, t/T\ < 1. 

The observations of the power-law scaling and of the 
random phases above / > f co validate the comparison of 
the continuum spectral range with the WTT. We test the 
scaling of the spectral energy with the energy flux P given 
by Eq. ([3]). In our experiments the energy input is pro- 
portional to the square of the wave amplitude, which is 
proportional to the square of the acceleration above the 
threshold of the parametric instability, P cx Aa 2 . Fig- 
ure[4]ja) shows the total energy in the spectrum as a func- 
tion of Aa 2 . The spectral power integrated from to 2 kHz 
varies as \[P. This integral is however dominated by the 
exponential and by the ZFSB ranges. To test the scaling 
of the spectral energy in the power-law range we compute 
JotkHz E{f)df. Fig. [^b) shows that this integral scales 
linearly with Aa 2 , contradicting the E u ~ P 1 / 2 prediction 
of the WTT. A similar observation of the linear Ej vs P 
dependence was already reported in experiments on sur- 
face waves by Falcon et al. [7] , which cast a doubt on the 
agreement between experiments and the WTT. In that ex- 
periment waves were excited at rather low frequency and 
the spectra showed the power-law scalings over the entire 
spectral range. 

Finally, to illustrate that our results do not depend on 
the container size, we performed experiments at much 
higher excitation frequency, /o =3.2 kHz. In this case, 
a circular container of the diameter d =180 mm accom- 
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Fig. 4: Spectral energy of capillary waves contained in differ- 
ent frequency ranges versus the energy flux injected by forcing 
(P oc Aa 2 ): (a) for the entire spectral range from to 2kHz 
and (b) in the power-law range, from 500Hz to 2kHz. 



modates 160 wave lengths on the radius (first harmonic, 
fi =1.6 kHz, A ~ 0.56 mm), substantially more than at 
/i =30 Hz shown in Fig. [T] The spectrum for the case of 
the high-frequency excitation is shown in Fig. [5j It has all 
the features observed at low frequency: discrete spectrally 
broadened peaks, exponential fit for the integrated spec- 
tral energy, and the ZFSB. Note that the spectral width 
of the ZFSB extends up to 200 Hz at high drive. Broad- 
ening of the harmonics is quite similar to that at low fre- 
quency, showing characteristic exponential tails (triangu- 
lar line shapes on the log-linear scale), Fig.[5jb). Thus, we 
conclude that the finite-box size effects are not important 
in these experiments above the threshold of modulation 
instability. 

Summarizing, we have shown that modulation instabil- 
ity plays a major role in the transition from monochro- 
matic parametrically excited wave to the broadband spec- 
tra. The instability is responsible for the strong mod- 
ulation of the wave amplitude and the breaking of the 
continuous waves into envelope solitons. This leads to a 
substantial shortening of the autocorrelation time and the 
correlation length which results in the detachment of the 
wave field from the container walls. At the same time, this 
very process leads to the spectral broadening of the wave 
harmonics which start overlapping to generate the wave 
continuum. In the continuum three-wave interactions be- 
tween two high-frequency waves and one low-frequency 
(/i) harmonic become possible. Shorter correlation times 
of the wave harmonics due to the MI randomize phases in 
the wave-wave interactions, leading to the conditions as- 
sumed in the WTT. When these conditions are satisfied, 
a power-law continuum is observed. 

Three distinct ranges in the spectra of the parametri- 
cally driven capillary waves have been identified. Spectral 
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Fig. 5: (a) Spectrum in the regime of high-frequency excitation 
at fo = 3.2 kHz, high forcing at Aa = 65g. A constant in 
the exponential fit A =0.7. (b) Zoomed-in part of the same 
spectrum. Dashed red spectral line corresponds to low forcing 
at Aa = 3g. 

energy contained in multiple harmonics obeys the expo- 
nential fit Ef = Eie~ A f/f 1 with A slowly decreasing with 
the increase in forcing from A = 1.2 to A = 0.8. Wave 
harmonics spectrally broaden as the forcing is increased. 
This broadening is observed for the first subharmonic and 
for its harmonics very close to the threshold of paramet- 
ric instability. The broadening increases proportionally to 
the forcing (acceleration), while the spectral line shapes 
remain unchanged, E n oc sech[b n (f — /„)]. 

Above some critical broadening of the wave harmonics, 
a continuum develops in the spectrum, starting at higher 
frequencies. This continuum is characterized by a power- 
law scaling, which at the highest drive shows Et ~ /~ 2 ' 8 , 
close to the one predicted by the WTT. However, the pre- 
diction of the dependence of the spectral energy on the in- 
jected energy flux differs from the theoretically predicted 
Ef cx P 1 / 2 : in experiments Ef scales linearly with the in- 
jected power P. Such a discrepancy may be due to the fact 
that the energy delivered into the power-law continuum is 
dissipated there, rather than being cascaded further down 
the spectrum. In other words, the assumption of the iner- 
tial interval in the continuum is not satisfied. 

We have found a criterion for the transition from the 
exponential scaling of Ef to the power law scaling. This 
occurs at the frequency of the wave harmonics whose av- 
erage length (time) of the modulation envelope becomes 
shorter than the period of the first harmonic, at r n <T\. 
This criterion corresponds to a strong reduction in the 
nonlinear interaction time between the first and the n-th 
wave harmonics, or in other words, leads to the phase- 
randomization in wave- wave interactions. 

Results obtained in this paper are not dependent on the 



size of the container: this has been checked by varying the 
frequency (and wave length A) of the first harmonic by a 
factor of more than 50 (16 for A) at the given size of the 
container. 

* * * 
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